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ABSTRACT 



We report an analysis of the interstellar 7-ray emission from the Chamaeleon, 
R Coronae Australis (R CrA), and Cepheus and Polaris flare regions with the Fermi 
Large Area Telescope. They are among the nearest molecular cloud complexes, within ~ 
300 pc from the solar system. The 7-ray emission produced by interactions of cosmic- 
rays (CRs) and interstellar gas in those molecular clouds is useful to study the CR 
densities and distributions of molecular gas close to the solar system. The obtained 
7-ray emissivities above 250 MeV are (5.9 ± 0.1 sta t iiosys) x 10~ 27 photons s _1 sr _1 
H-atom" 1 , (10.2 ± 0.4 stat t[^ By8 ) x 10~ 27 photons s" 1 sr^ 1 H-atom^ 1 , and (9.1 ± 
0.3 s tat lo6sys) x 10~ 27 photons s^ 1 sr _1 H-atom -1 for the Chamaeleon, R CrA, and 
Cepheus and Polaris flare regions, respectively. Whereas the energy dependences of the 
emissivities agree well with that predicted from direct CR observations at the Earth, 
the measured emissivities from 250 MeV to 10 GeV indicate a variation of the CR den- 
sity by ~ 20 % in the neighborhood of the solar system, even if we consider systematic 
uncertainties. The molecular mass calibrating ratio, Xqq = A r (H 2 )/T4 / co7 is found to 
be (0.96 ± 0.06 stat taiisys) xl 20 H 2 -molecule cm" 2 (K km s -1 ) -1 , (0.99 ± 0.08 sta t 
loaosys) xl ° 20 H 2 -molecule cm" 2 (K km s" 1 )" 1 , and (0.63 ± 0.02 stat t^'^eys) xlO 20 
H2-molecule cm -2 (K km s^ 1 ) -1 for the Chamaeleon, R CrA, and Cepheus and Polaris 
flare regions, respectively, suggesting a variation of Aco in the vicinity of the solar 
system. From the obtained values of AcO; the masses of molecular gas traced by Wco 
in the Chamaeleon, R CrA, and Cepheus and Polaris flare regions are estimated to be 
~ 5xl0 3 Mq, ~ 10 3 Mq, and ~ 3.3 xlO 4 M©, respectively. A comparable amount of 
gas not traced well by standard H 1 and CO surveys is found in the regions investigated. 
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1. Introduction 



Observations of high-energy 7-ray emission (E > 30 MeV) from molecular clouds can be used 
to study the cosmic-ray (CR) production, the CR density, and the distribution of the interstel- 
lar medium (ISM) in such systems. 7-rays are produced in the ISM by interactions of high- 
energy CR protons and electrons with the interstellar gas, via nucleon-nucleon collisions, electron 
Bremsstrahlung, and inverse Compton (IC) scattering. Since the 7-ray production cross section is 
almost independent of the chemical or thermodynamic state of the ISM, and the interstellar gas is 
essentially transparent to those high-energy photons, observations in 7-rays have been recognized as 
a powerful probe of the distribution of interstellar matter. If the gas column densities are estimated 
with good accuracy by observations in other wavebands such as radio, infrared, and optical, the 
CR spectrum and density distributions can be examined as well. Molecular clouds that are within 
1 kpc from the solar system (namely nearby molecular clouds) and have masses greater than a few 
10 3 Mq are well suited for an analysis of their 7-ray emission to investigate the distribution of CR 
densities and interstellar gas since they are observed at high latitudes and therefore largely free 
from confusion with the strong emission from the Galactic plane. Study of such nearby molecular 
clouds in 7-rays can be dated back to the COS-B era (e.g., Bloemen et al. 1984) and was advanced 
by the EGRET on board Compton Gamma-Ray Observatory (e.g., Hunter et al. 1994). Although 
some important information has been obtained on properties of CRs and the ISM by these early 
observations, detailed studies have only been performed on giant molecular clouds with masses 
greater than ~ 10 5 Mq such as the Orion complex (e.g., Digel et al. 1999). The data above 1 GeV, 
which are crucial to study CR nuclei spectra, suffered from the limited photon statistics, angular 
resolution, and energy coverage of these early missions. 

The advent of the Fermi Gamma-ray Space Telescope launched in 2008 has improved the 
situation significantly. The sensitivity of the LAT (Large Area Telescope) on board Fermi is more 
than an order of magnitude better than that of the EGRET, and enables resolving more point 
sources and studying the diffuse 7-ray emission with unprecedented sensitivity. In addition, newer 
surveys of the ISM (e.g., Dame et al. 2001, Kalberla et al. 2005, and Grenier et al. 2005) allow us 
to investigate the CR spectral and density distributions with better accuracy. 

Here we report a Fermi LAT study of diffuse 7-rays from the Chamaeleon, R Coronae Australis 
(R CrA), and Cepheus and Polaris flare molecular clouds. They are among the nearest (< 300 pc 
from the solar system) molecular clouds exhibiting star formation activity. Although EGRET 
observed 7-ray emission associated with the molecular gas in the Chamaeleon region (Grenier et 
al. 2005), no detailed study of CR and matter distributions for the Chamaeleon and R CrA regions 
has been performed yet since they have rather small masses (< 10 4 M Q , about 1/10 of that of the 
Orion molecular cloud) and consequently small 7-ray fluxes. We also analyzed in detail the region 
of the Cepheus and Polaris flares which was included in the Fermi LAT study of the second Galactic 
quadrant (Abdo et al. 2010b). It is located in the direction almost opposite to the Chamaeleon 
region in the Gould Belt (see, e.g., Perrot and Grenier 2003), therefore we can investigate the 
distribution of the CR density over several hundred pc but still inside the coherent environment 



of the Gould Belt. This paper is organized as follows. We first describe the observations as well 
as the sky model preparation and the data analysis in Section [21 and show the obtained results in 
Section [3l We then discuss the CR and matter distributions in Section [5] and give conclusions in 
Section [5j 

2. Data Analysis 

2.1. Observations and Data Reduction 

The LAT, on board the Fermi Gamma-ray Space Telescope, is a pair-tracking detector to 
study 7-rays from ~ 20 MeV to more than 300 GeV. It consists of an array of 4x4 conversion 
and tracking modules built with tungsten foils and silicon microstrip detectors to measure the 
arrival directions of incoming 7-rays and a hodoscopic cesium iodide calorimeter to determine the 
photon energies. The modules are surrounded by 89 segmented plastic scintillators serving as an 
anticoincidence detector to reject charged-particle background events. A detailed description of the 
LAT instrumentation can be found in Atwood et al. (2009) and the on-orbit calibration is discussed 
in Abdo et al. (2009a). 

Science operations with the LAT started on 2008 August 4. For this analysis we have accu- 
mulated events obtained from 2008 August 4 to 2010 May 9. During this time interval the LAT 
was operated in sky survey mode nearly all of the time and scanned the 7-ray sky with relatively 
uniform exposure over time (within 10% in regions studied). We used the standard LAT analysis 
software, Science ToolaS version v9rl6p0 and the response function P6_V3_DIFFUSE, which was 
developed to account for the detection inefficiencies due to pile-up and accidental coincidence of 
events (Rando et al. 2009). We applied the following event selection criteria to the 7-ray events: 
(1) events must satisfy the standard low-background event selection (so-called diffuse class events; 
Atwood et al. 2009), (2) the reconstructed zenith angles of the arrival direction of photons are 
less than 100° in order to reduce contamination of photons from the bright Earth rim, and (3) 
the center of the LAT field of view is within 52° from the zenith direction of the sky, in order to 
exclude the data obtained during the relatively short time intervals of pointed observations when 
the rocking angle of the LAT was larger. The exposure maps were generated with the same 100° 
limit on zenith angle for each direction in the sky. 

The count maps obtained (E > 250 MeV) in the Chamaeleon, R CrA, and Cepheus and 
Polaris flare regions are shown in Figure [TJ we set the lower energy limit at 250 MeV to utilize 
good angular resolution (68% containment radius is < 1.5° above 250 MeV) and the upper energy 
limit at 10 GeV because of limited photon statistics. We also show positions of point sources 



Available from the Fermi Science Support Center ( http : / /fermi .gsfc . nas a.gov /ssc/ 1 . 



with high significance (test statistic, T£q greater than 50) and 2.6 mm carbon-monoxide CO line 
intensities on the maps. The yellow square indicates the region of interest (ROI) analyzed for each 
of the regions. In order to take into account the spillover from point sources outside of the ROIs, 
we also included point sources lying just outside (< 5°) of the region boundaries. Contamination 
due to the diffuse emission from the interstellar gas outside of the ROIs is also taken into account 
through the convolution of maps larger than the ROIs. 

2.2. Model Preparation 

Since the ISM is optically thin to 7-rays in the energy range considered in the paper and the 
CR spectrum is not expected to vary significantly within small regions, the 7-ray intensity from 
CR protons and electrons interacting with the interstellar gas may be modeled as a sum of emission 
from separate gas phases (e.g., Lebrun et al. 1983). This approach has been successfully applied in 
recent studies of diffuse 7-rays by the LAT (e.g., Abdo et al. 2010b and Ackermann et al. 2011). 
We followed this method and prepared template maps as described below. 

2.2.1. H 1 and CO Maps 

We calculated the column densities iV(H 1) of atomic hydrogen from the Leiden/ Argentine/Bonn 
Galactic H 1 survey by Kalberla et al. (2005). The optical depth correction of the H 1 gas is applied 
under the assumption of a uniform spin temperature To, = 125 K, the value which has often been 
used in previous studies (e.g., Abdo et al. 2009c and Abdo et al. 2010b). This choice of T$ allows 
us to directly compare our results with other studies. In order to evaluate systematic uncertainties 
due to the optical depth correction, we also tried several different choices of Ts as described in 
Section [3l We note that the true Ts is likely not to be uniform even in small regions like the ones 
we are studying, but exploring a non-uniform Ts is beyond the scope of this paper. We separated 
the Hi column densities in two regions along the line of sight; one corresponds to the local region 
(< 300 pc) to take into account the ambient atomic gas surrounding the molecular cloud and the 
other corresponds to the rest to take into account the remaining gas along the line of sight. From 
the velocity distribution of the CO emission which traces the molecular cloud, we determined the 
boundary as shown in Figure [2j The local velocity range is —10 < vlsr < 15 km s _1 for the 
Chamaeleon region and — 15 < vlsr < 15 km s" 1 for the R CrA region for b > —10°. Below —10°, 
I^lsrI is increased to 80 km s _1 at b = —20.5° in both regions, since the Hi gas at such high 
latitude is likely to be local. The obtained iV(H 1) maps are shown in the top panels of Figure [3] 
and Figure U] for the Chamaeleon region and the R CrA region, respectively. For the Cepheus and 
Polaris flare region, the cut falls in between the Gould Belt lines and the Local-Arm lines at —8 



2 TS is defined as TS = 2(lnL — lnLo), where L and Lo are the maximum likelihoods obtained with and without 
the source included in the model fitting, respectively; see Mattox et al. (1996). 
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Fig. 1. — 7-ray count maps above 250 MeV for the Chamaeleon (top left), R CrA (top right), and 
Cepheus and Polaris flare (bottom left) regions, smoothed with a Gaussian of a standard deviation 
a = 0.5° for display. The contours indicate intensities Wqo of the 2.6 mm line of CO (with the 
levels of 4, 8, 12, and 16 K km s" 1 ) by Dame et al (2001), as a standard tracer of the molecular 
gas. Cyan circles show the positions of point sources with high significance (TS > 50) in the First 
Fermi LAT catalog (1FGL) by Abdo et al. (2010a). The yellow squares indicate the ROI analyzed 
for each of the regions. Point sources outside of this ROI but inside the cyan square are taken into 
account in the analysis. 



km s at b < 15°, and then opens up to —100 km s at b = 24°. Since the amount of Hi gas in 
the Local-Arm and beyond is comparable to that in the Gould Belt for the Cepheus and Polaris 
flare region in 15° < b < 20° (see also Figure [5]), using two Hi template maps is crucial. For the 
Chamaeleon and the R CrA regions, the non-local H I gas is almost negligible. 

The integrated intensities of the 2.6 mm line, Wco, have been derived from the composite 
survey of Dame, Hartmann, & Thaddeus (2001). We used this Wco ma P as a standard molecular- 
gas tracer. For better signal-to-noise ratio, the data have been filtered with the moment-masking 
technique (Dame 2011) to reduce the noise while keeping the resolution of the original data. Since 
most of the molecular gas turned out to be local according to our velocity cuts, we used only the 
local CO map in the 7-ray analysis. 

2.2.2. Excess Av Map 

Dust is a commonly-used tracer of the neutral interstellar gas. By comparing the 7-ray ob- 
servations by EGRET with radio surveys and the dust thermal emission, Grenier et al. (2005) 
reported a considerable amount of gas at the interface between the atomic/molecular phases in 
the solar neighborhood, associated with cold dust but not properly traced by H 1 and CO surveys. 
This finding was confirmed by dedicated analyses of the diffuse 7-ray emission with the Fermi LAT 
(Abdo et al. 2010b and Ackermann et al. 2011). In order to take into account this additional 
interstellar gas, we constructed visual extinction (Av) maps, based on the extinction maps derived 
by Schlegel et al. (1998) from IRAS and COBE DIRBE data. The Av map on the assumption of a 
constant gas-to-dust ratio provides an estimate of the total column densities. After fitting a linear 
combination of the iV(H 1) and Wco maps through a minimum sum-of-square-residuals criterion in 
each ROI separately, we thus obtained a residual extinction map, Av rcs map, accounting for the 
additional gas which is not properly traced by the Hi and CO surveys. Negative residuals are likely 
to be due to the fluctuation of the original Av map. For simplicity, we clipped data around in 
the Av res map. 

We present gas maps used for the analysis of the Chamaeleon, R CrA, and Cepheus and Polaris 
flare regions in Figures El [H and El respectively. 
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Fig. 2. — Latitude-velocity diagrams of the intensity of the 21 cm line (log scale in atom cm -2 ) for 
the Chamaeleon (top left), R CrA (top right), and Cepheus and Polaris flare (bottom left) regions. 
The dashed lines indicate the region boundary between the local region (< 300 pc) and the rest. 
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Fig. 3. — Template gas maps for the Chamaeleon region; iV(Hi) (local) (top left) and iV(Hi) (non- 
local) (top right) in units of 10 20 H-atoms cm" 2 , Wco (bottom left) in units of K km s _1 , and 
Av res (bottom right) in units of magnitudes. The two iV(Hi) maps have been smoothed with a 
Gaussian of a standard deviation a = 1° for display while the other two maps have been smoothed 
with a Gaussian of a = 0.25° in order to keep fine structures seen in Wqq and Av res distributions. 
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Fig. 4. — Template gas maps for the R CrA region; N(H i) (local) (top left) in units of 10 20 atoms 
cm -2 , Wqo (bottom left) in units of K km s _1 , and ^4v rcs (bottom right) in units of magnitudes. 
They have been smoothed in the same way as maps in Figure El The non-local iV(H i) is almost 
in our ROI and hence is not shown here. 





1lS 114 !Si 

Gnlarr*: LiniglEuita Vim) 




1lS 11 » 



Fig. 5. — Template gas maps for the Cepheus and Polaris flare region; N(H i) (local) (top left) and 
./V(Hl) (non-local) (top right) in units of 10 20 atoms cm" 2 , Wco (bottom left) in units of K km 
s _1 , and Av rcs (bottom right) in units of magnitudes. They have been smoothed in the same way 
as maps in Figure 



2.2.3. IC, Point Sources and Isotropic Component 



In addition to the gas-related components described above, we also need to take into account 
7-rays from IC scattering, contributions of point sources, extragalactic diffuse emission, and in- 
strumental residual background. To model the IC emission, we used GALPROP (e.g., Strong & 
Moskalenko 1998), a numerical code which solves the CR transport equation within our Galaxy and 
predicts the 7-ray emission produced via interactions of CRs with the ISM. IC emission is calculated 
from the distribution of propagated electrons and the radiation field model developed by Porter et 
al. (2008). Here we adopt the IC model map produced in the GALPROP run 54_77Xvarh7^| as 
a baseline model, in which the CR electron spectrum is adjusted based on the Fermi LAT mea- 
surement. In this model, the CR source distribution model is adjusted to the LAT data, and is 
somewhat more concentrated to the inner Galaxy than the pulsar distribution by Lorimer (2004). 
To take into account uncertainties of the CR electron spectrum and radiation field on the Galactic 
scale, we set the normalization of this IC component free in each energy bin when we perform the 
fit (see Section T2.3p . In order to evaluate the systematic uncertainty due to IC models, we also tried 
four other models, which are constructed under different assumptions about the distribution of CR 
sources such as supernova remnants (Case & Bhattacharya 1998) and pulsars (Lorimer 2004), and 
intensity of the interstellar radiation field depending on the input luminosity of the Galactic bulge 
component (e.g., Ackermann et al. 2012). 

To take into account 7-ray point sources, we used the 1FGL catalog based on the first 11 
months of the science phase of the mission (Abdo et al. 2010a). We included in our analysis point 
sources detected with TS > 50 in the 1FGL catalog and other significant point sources included in 
the 2FGL catalog as described in Section 12.31 

To represent the sum of the extragalactic diffuse emission and the residual background from 
the misclassified CR interactions in the LAT, we adopted a publicly available isotropic spectrum^ 
obtained by a fit to emission from the high latitude sky (b > 30°). This component is fixed in our 
analysis. The uncertainty due to this isotropic term will be discussed in Section [3l 



http://www.mpe.mpg.de/~aws/propagate.html 



4 isotropic_iem_v02.txt from |http: / / fermi.gsfc.nasa.gov / ssc / data/access /lat /BackgroundModels.html 



2.3. Analysis Procedure 



With the usual assumptions of optical thinness and that CRs uniformly thread the ISM, 7-ray 
intensity I~,(l,b) (s _1 cm -2 sr _1 MeV -1 ) at a given energy can be modeled as 

2 

I 7 {l,b) = ^2qm,i-NQli)(l,b)i + qco-Wco(l,b) 

i=l 

+ q Av -Av TCS (l,b)+c lc -hc(l,b)+I is o + ^2~PSj (1) 

3 

where sum over i represents the two regions (local and non-local regions), qua (s -1 sr" 1 MeV -1 ), 
Qco (s" 1 cm" 2 sr" 1 MeV" 1 (K km s" 1 )" 1 ), and gAv (s" 1 cm" 2 sr" 1 MeV" 1 mag" 1 ) are the 
emissivity per Hi atom, per Wco unit, and per Av Tes magnitude, respectively. lie an d ii so are 
the IC model and isotropic background intensities (s" 1 cm" 2 sr" 1 MeV -1 ), respectively, and PSj 
represents contributions of point sources. Hard 7-ray emission with a characteristic bubble shape 
above and below the Galactic center (usually called the "Fermi bubbles ") that was found in the 
Fermi LAT data (e.g., Su et al. 2010) has large spatial extent in the R CrA region. As described 
below, we found residuals with a hard 7-ray spectrum in the R CrA region, and included an 
additional template in equation (TT]). 

The 7-ray data in our ROIs were binned in 0.25° x 0.25° pixels and fitted with equation ([1]) in 
8 logarithmically equally-spaced energy bins from 250 MeV to 10 GeV using a binned maximum- 
likelihood method with Poisson statistics. Low photon statistics and poor angular resolution at 
low energy (~ 1.5° at 250 MeV under the 68% containment radius) do not allow us to separate 
components reliably. For convolution of diffuse emission with the instrumental response functions, 
we assumed an E~ 2 spectrum and the integrated intensities were allowed to vary in each of the 
8 energy bins. Changing the fixed spectral shape index over the range from —1.5 to —3.0 has 
negligible effect on the obtained spectrum. Data in the R CrA region with energies above 4 GeV 
are grouped in a single bin to get larger statistics. 

We started the analysis for the Chamaeleon region with point sources detected with high 
significance (TS > 100) in the 1FGL catalog. The normalizations for each energy bin are allowed to 
vary for sources inside our ROI. We also included sources lying outside (< 5°) ROI, with the spectral 
parameters fixed to those in the 1FGL catalog. We first fitted the model to LAT data without the 
Av res map, and then included it and confirmed that the likelihood improved significantly; the test 
statistic, defined as TS = 2(lnLi — InLo), where Lq{L\) is the likelihood without(with) additional 
component, is 704 for 8 more free parameters for the energy range from 250 MeV to 10 GeV. 
Figure [6] shows residual (data minus fitted model) map obtained from the fit without the ^4v res 
map. Positive residual counts are seen where we have positive Av res (Figure [3]). We thus confirmed 
that the positive Av rcs traces the gas not well measured by H 1 and CO surveys, and included 
the Av res map in the following analysis. We repeated the same procedure and obtained the same 
conclusion for the Cepheus and Polaris flare region. In the R CrA region, large residual clumps are 



still seen around (0° < I < 15°) and (—30° < b < —15°) even if we included the ^4v res map as shown 
by Figure probably due to the southern Fermi bubble. In order to accout for these residuals, 
we used a flat template model map with the shape as shown in Figure [7J with a free normalization 
in each energy bin. We note that the template map is just to accommodate the residuals not to 
investigate the Fermi bubble. The residuals are improved significantly as shown in the bottom 
right panel of Figure [H and the intensity of this template is too low to significantly impact on the 
local HI gas emissivity as shown in the top right panel of Figure \12[ 

We then lowered the threshold for point sources down to TS = 50. Although the fit improves in 
terms of the log-likelihood, the effect on the emissivities associated with gas maps is negligible for the 
three regions (smaller than the statistical error). However, some point-like excesses corresponding to 
objects not included in the 1FGL catalog are seen in the R CrA region (2FGL J1830. 2-4441, 2FGL 
J1816.7-4942, and 2FGL J1825. 1-5231) and Cepheus and Polaris flare region (2FGL J2022. 5+7614 
and 2FGL J2009.7+7225). These may be sources that became luminous after the 1FGL catalog 
was published. We thus included those sources and confirmed that gas emissivities are almost 
unaffected while the residual map becomes flatter. We therefore adopted the model described 
by equation ([1]) (plus a flat template for the R CrA region) with point sources with TS > 50 in 
the 1FGL catalog (plus additional sources in the R CrA, and Cepheus and Polaris flare regions 
described above) as our baseline model. We note that an unassociated source in the 1FGL catalog, 
1FGL J1903.8-3718c, is located on a CO core of the R CrA molecular cloud and was recognized in 
the 1FGL catalog as potentially being spurious. We thus did not include this source in our model. 
We confirmed that the obtained emissivities are almost the same if we mask the region of the CO 
core. 
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Fig. 6. — 7-ray residual map obtained from the fit without the Av rcs map in units of standard 
deviations above 250 MeV for the Chamaeleon region smoothed with a Gaussian of a = 0.5°. 
Positive residual counts are seen where we have positive Av res (Figure [3]). 




Fig. 7. — 7-ray residual map obtained from the fit in units of standard deviations above 250 MeV for 
the R CrA region smoothed with a Gaussian of a = 0.5°. Dashed lines indicate region boundaries 
for the additional flat template map. 



3. Results 



Figures [HJ [H and [10] show the 7-ray data count maps, fitted model count maps and the residual 
maps for the Chamaeleon, R CrA, and Cepheus and Polaris flare regions, respectively, in which 
the residuals are expressed in units of approximate standard deviations (square root of the model 
counts). The residual maps show no conspicuous structures, indicating that our model reasonably 
reproduces the data, particularly the diffuse emission. For illustrative purposes, we present the 
fitted model count maps for the Chamaeleon region decomposed into each gas component in Figure 
[TT| in which the 7-ray emission from Hi, that from the molecular gas traced by Wqo, an d that 
inferred from the ^4v res map are shown. Although the distribution of iV(Hl) is rather uniform in 
our ROI, it exhibits some structures and allows us to derive the emissivity of the H 1 gas, which 
is proportional to the flux of ambient CRs. The distribution of Wco is highly structured and is 
concentrated in the longitude range from 295° to 305° and in the latitude range from —12° to —20°. 
The gas traced by ^4v res lies at the interface between the H 1 component (atomic gas) and the CO 
component (molecular gas), and has a mass (proportional to the 7-ray counts) comparable to or 
larger than that of the molecular gas traced by Wco- 

Figure [T2] shows the fitted spectra for each component. Although the contributions from IC 
and isotropic components are large, the spectra of each gas component are reliably constrained due 
to their characteristic spatial distributions as shown in Figures [3] [U and [5j The hard spectra of 
the IC term and the flat template model component of the R CrA region are likely to be due to 
the southern Fermi bubble. 

The fit results are summarized in Tables HJ [2j and [3] for the case of Ts = 125 K in the 
Chamaeleon, R CrA, and Cepheus and Polaris flare regions, respectively. The integrated H 1 emis- 
sivity for the Chamaeleon region above 250 MeV is (5.9 ± 0.1 stat t^osys) x 10" 27 photons s" 1 
sr" 1 H-atom" 1 , and those of the R CrA, and Cepheus and Polaris flare regions are (10.2 ± 0.4 stat 
-lifsys) x 10" 27 photons s" 1 sr" 1 H-atom" 1 and (9.1 ± 0.3 sta t to^sys) x 10 ~ 27 photons s" 1 sr" 1 
H-atom -1 , respectively. (See below for the evaluation of the systematic uncertainty from the sky 
model.) 

Figure [13] shows the emissivity spectra of each gas component in the Chamaeleon region under 
the assumption of Ts = 125 K. In order to examine the systematic uncertainty due to the optical 
depth correction, we also tried to fit the data with maps obtained by assuming Ts = 100 K 
and under the approximation that the gas is optically-thin. We evaluated the uncertainty of the 
isotropic component to be ±10% by comparing the model we adopted and those derived in other 
LAT studies of mid-latitude regions (Abdo et al. 2009b and Abdo et al. 2009c). We thus reran the 
analysis described in Section 12.31 assuming a 10% higher and lower intensity for the fixed isotropic 
component. We also investigated the effect on the systematic uncertainty due to the IC component 
by using different IC model maps, as described in Section 12.2.31 The effects of the uncertainty 
of the Ts, isotropic component, and IC models are quite comparable, therefore we added them. 
The obtained systematic uncertainty is comparable to or slightly larger than the statistical error 
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Fig. 8. — Data count map (top left), fitted model count map (top right), and residual map (bottom 
left) in units of standard deviations above 250 MeV under the assumption of Tg = 125 K for the 
Chamaeleon region. The residual map has been smoothed with a Gaussian of a = 0.5°. 
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Fig. 9. — The same as Figure [8] for the R Cr A region. 




Fig. 11. — Fitted model count maps above 250 MeV for the Chamaeleon region; Hi component 
(top left), CO component (top right) and Av res component (bottom left). 
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Fig. 12. — 7-ray spectra of each component (each gas phase, IC, isotropic component, point sources, 
and flat template model for the R CrA region) for the Chamaeleon (top left), R CrA (top right), 
and Cepheus and Polaris flare region (bottom left). The Hi (non-local) for the R CrA region and 
the IC component for the Cepheus and Polaris flare are not shown here since they are negligible. 



Table 1: Gas emissivities in the Chamaeleon regions with their statistical uncertainties (assuming 
Ts = 125 K for the Hi maps preparation). 





Out i 
y 111,1 




yAv 


(GeV) 








0.25-0.40 


2.3±0.1 


0.30±0.06 


0.57±0.03 


0.40-0.63 


1.50±0.06 


0.33±0.04 


0.34±0.02 


0.63-1.00 


0.99±0.04 


0.17±0.02 


0.21±0.01 


1.00-1.58 


0.61±0.03 


0.13±0.01 


0.125±0.008 


1.58-2.51 


0.29±0.02 


0.069±0.009 


0.060±0.006 


2.51-3.98 


0.13±0.01 


0.027±0.006 


0.033±0.004 


3.98-6.31 


0.08±0.01 


0.010±0.004 


0.014±0.003 


6.31-10.00 


0.019±0.006 


0.007±0.003 


0.005±0.002 


total 


5.9±0.1 


1.04±0.08 


1.36±0.04 


Notes, units; gHl,i(10" 


- 27 s- 1 sr" 1 ), gco(10" 6 


cm -2 s _1 sr -1 (K km s — 1 )" 


_1 ), gAv(10 — 5 cm -2 s —1 sr -1 mag -1 ) 



as shown by Figure [13] for qm (top left panel). On the other hand, the systematic uncertainty of 
the qco and gAv is smaller than the statistical error and is not shown in this figure for clarity. The 
obtained spectra in the R CrA region and the Cepheus and Polaris flare region are summarized in 
Figures ?? and [151 respectively. We performed the same procedure and evaluated the systematic 
uncertainties due to the Ts, isotropic component, and IC models as a shaded area in the figures. 

Evaluation of the systematic uncertainty is crucial for this study. We therefore performed two 
more tests. 

1. We modified the longitude range from -20° < I < 20° to -30° < I < 10° for the R CrA 
region, and found that the obtained Hi emissivity was lower by ~ 10 %. This is due to a 
coupling between the H I and IC components, and take this effect on the H I emissivity into 
account in the evaluation of the overall systematic uncertainty. 

2. The prominent Hi cloud at 280° < I < 295° and -28° < b < -20° is not spatially associated 
with the Chamaeleon molecular cloud, and it may not be physically related to the molecular 
cloud. We masked this region and found that the H I emissivity is almost unaffected. 

The resultant peak-to-peak uncertainty of the local H I emissivity is less than ~20 % across 
the energy range for three regions investigated. 



Table 2: The same as Tabled] for the R CrA region. 



Energy 


qm 


<?CO 


<7Av 


(GeV) 








0.25-0.40 


4.1±0.3 


0.7±0.1 


1.0±0.2 


0.40-0.63 


2.6±0.2 


0.45±0.08 


0.6±0.1 


0.63-1.00 


1.6±0.1 


0.37±0.05 


0.41±0.06 


1.00-1.58 


0.96±0.09 


0.23±0.03 


0.15±0.04 


1.58-2.51 


0.53±0.06 


0.10±0.02 


0.10±0.03 


2.51-3.98 


0.25±0.04 


0.05±0.01 


0.06±0.02 


3.98-10.00 


0.14±0.03 


0.019±0.008 


0.008±0.001 


total 


10.2±0.4 


1.9±0.2 


2.3±0.2 


Notes, units; g H i(10~ 27 


s" 1 sr" 1 ), «j C o(10- 6 


cm -2 s" 1 sr" 1 (K km s^ 1 )" 1 ), g A v(10 -5 


cm~ 2 s" 1 sr -1 mag -1 ) 



Table 3: The same as Table Q] for the Cepheus and Polaris flare region. 



Energy 
(GeV) 


<?HI 


qco 


<?Av 


0.25-0.40 


3.5±0.1 


0.52±0.04 


0.53±0.06 


0.40-0.63 


2.37±0.06 


0.29±0.02 


0.41±0.04 


0.63-1.00 


1.6±0.3 


0.21±0.01 


0.22±0.02 


1.00-1.58 


0.97±0.03 


0.105±0.008 


0.11±0.02 


1.58-2.51 


0.48±0.02 


0.061±0.005 


0.06±0.01 


2.51-3.98 


0.20±0.01 


0.029±0.003 


0.030±0.006 


3.98-6.31 


0.11±0.01 


0.012±0.002 


0.016±0.004 


6.31-10.00 


0.035±0.006 


0.006±0.001 


0.006±0.003 


total 


9.2±0.3 


1.23±0.05 


1.38±0.08 



Notes, units; <jhi(10 27 s 1 sr 1 ), <jco(10 6 cm 2 s 1 sr 1 (K km s 1 ) 1 ), (jAv(10 5 cm 2 s 1 sr 1 mag 1 ) 
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Fig. 13. — Emissivity spectrum of the local Hi gas (top left), that per Wco unit (top right) and 
that per unit Av rcs magnitude (bottom left) of the Chamaeleon region. The shaded area shows 
systematic uncertainties for Hi (see text for details). 




Fig. 14. — The same as Figure [13] for the R Cr A region. 
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Fig. 15. — The same as Figure [13] for the Cepheus and Polaris flare region. 



4. Discussion 



4.1. CR Density and Spectrum Close to the Solar System 

It is generally believed that supernova remnants are the primary sources of the Galactic CRs. 
Due to strong selection effects in detecting supernova remnants, their spatial distribution is not well 
determined. Therefore a smooth, axisymmetric distribution has often been assumed in theoretical 
calculations, resulting in a smooth decline of CR density as a function of the Galactocentric radius. 
This assumption, however, should be tested against observations. In Figure [161 we compare the 
obtained H I emissivity spectra in the Chamaeleon, R CrA, and Cepheus and Polaris flare regions 
(for Ts = 125 K) and model emissivity spectra^ for the local interstellar spectrum (LIS) used in 
Abdo et al. (2009c), based on local direct measurement of CRs. The spectral shapes for the three 
regions studied here agree well with the LIS models, indicating that the CR nuclei have similar 
spectral distribution in the vicinity of the solar system. On the other hand, the absolute emissivities 
differ among regions. The emissivities of the three regions studied differ by ~ 50%, estimated from 
the total qui shown in Tables Q3 [21 and [3l We note that the systematic uncertainty of the LAT 
effective area (5% at 100 MeV and 20% at 10 GeV; Rando et al. 2009) does not affect the relative 
value of emissivities among these regions in Figure [16j Although the emissivities of the R CrA 
region and the Cepheus and Polaris flare region are comparable, that of the Chamaeleon region 
is lower by ~ 20%, even if we take the systematic uncertainty into account (see Section [3]). As 
a further test, we fixed the emissivity of the Chamaeleon region to that of the model for the LIS 
with the nuclear enhancement factor of 1.84 and performed the fitting. The fit turns out to be 
significantly worse: the obtained ln(L) is lowered by 153 with 8 less free parameters. In addition, 
the normalization of the IC term is lowered by more than a factor of three, although we cannot 
rule out such a low IC flux (low CR electron flux) for the direction of the Chamaeleon region. The 
effect of unresolved point sources is small, since we have verified that the obtained emissivities are 
robust against the lower threshold for point sources between TS = 50 and 100 (see Section [2.3p . 
We also confirmed that the residual excess of photons around (/ = 280° to 288°, b = —20° to —12°; 
see the bottom panel of Figure [8]) does not affect the local H I emissivity very much. Thus the 
total systematic uncertainty of the Chamaeleon region is conservatively estimated to be ~ 15% at 
most (mainly due to the Tg, isotropic component and IC models), indicating a difference of the 
CR density between the Chamaeleon and the others as shown in Figure [161 

If the CR density has a variation by a factor of 1.2-1.5 in the neighborhood of the solar system, 
this requires a serious reconsideration of a smooth CR density often adopted for simplicity, and 
may have an impact on the study of the CR source distribution and diffuse 7-ray emission. We 
note that CR sources are stochastically distributed in space and time, and this may produce a CR 



The model is calculated from the LIS compatible with the CR proton spectrum measured by Alcaraz et al. (2000) 
and Sanuki et al. (2000), under the assumption that the nuclear enhancement factors (the correction terms to take 
into account the contribution from nuclei heavier than protons in both CRs and interstellar matter) are 1.45 and 1.84 
(Mori 2009); see Abdo et al. (2009c). 



anisotropy depending on the propagation conditions as discussed by, e.g., Blasi & Amato (2011a) 
and Blasi & Amato (2011b). Study of other regions and more detailed theoretical calculations will 
be needed to further investigate this issue. 



4.2. 



Molecular Masses in the Interstellar Clouds Studied 



Since the 7-ray production is almost independent of the chemical or thermodynamical state of 
the interstellar gas, the 7-ray observation is a powerful probe to investigate the molecular mass cali- 
bration ratio, Aco, defined as A^B^/Wco- Under the hypothesis that the same CR flux penetrates 
the Hi and CO phases of an interstellar complex, we can calculate Aco as Aco = Qco I (2<Zhi) , as 
shown in Figure [T71 (left). The linear relation supports the assumption that Galactic CRs penetrate 
these molecular clouds uniformly to their cores. This also indicates that any contamination from 
point sources and CR spectral variations in molecular clouds analyzed here is small. 

We have obtained Aco by fitting the relation between qm and qco with a linear function 
using a maximum-likelihood method; the Aco values are (0.96 ± 0.06 sta t loilLys) x ^® 2 ° cm_2 



(K km s- 1 )- 1 , (0.99 ± 0.08 stat Ig^sys) x 10 20 cm" 2 (K km a" 1 )" 1 , and (0.63 ± 0.02 stat ±°;^ sys ) 



x 10 20 cm" 2 (K km s" 1 )" 1 for the Chamaeleon, R CrA, and Cepheus and Polaris flare regions, 
respectively. The obtained value of Aco f° r the Cepheus and Polaris flare region is ~ 20 % lower 
than that reported by Abdo et al. (2010b). Abdo et al. (2010b) includes in their study also 
the Cassiopeia molecular cloud in the Gould Belt, and due to different ROIs considered, the 
emissivity was also different. Aco of the Chamaeleon region is similar to that of the R CrA region, 
whereas that of the Cepheus and Polaris flare region is ~ 2/3 of the others. The LAT data thus 
suggest a variation of Aco on a ~ 300 pc scale. 

We can estimate the CO-bright molecular mass for the Chamaeleon, R CrA, and Cepheus and 
Polaris flare regions. The mass of the gas traced by Wco is expressed as 



where d is the distance to the cloud, toh is the mass of the hydrogen atom and fx = 1.36 is the 
mean atomic weight per H-atom (Allen 1973). From this equation the mass of gas traced by CO 
is expressed as M C o in Tabled we obtained ~ 5xl0 3 M , ~ 10 3 M , and ~ 3.3xl0 4 M for the 
Chamaeleon, R CrA, and Cepheus and Polaris flare regions, respectively. The obtained mass of the 
Cepheus and Polaris flare region is ~ 20 % lower than that reported by Abdo et al. (2010b) due to 
the different value of Aco- Our estimates for the Chamaeleon and the R CrA regions are ~ 1/2 of 
those by Dame et al. (1987); they obtained ~10 4 M and ~ 3xl0 3 M and for the Chamaeleon 
and R CrA regions, respectively, under the assumption of a high value of Aco = 2.7 x 10 20 cm" 2 
(K km s -1 ) -1 . 

Using the relation between qhi and q^ v , we can also calculate the mass of additional interstellar 
gas traced by Av rcs . Figure [T71 shows the results of the fitting by a linear relation, q^ v = Aa v • Qm- 




(1) 



Chamaeleon region 
' R CrA region 
Cepheus and Polaris flare region 
model for LIS (e M = 1.84) 
model for LIS (e M = 1.45) 



10" 



10 J 

Energy (MeV) 



Fig. 16. — H I emissivity spectra of the Chamaeleon, R CrA, and Cepheus and Polaris flare regions 
compared with the model for the LIS with nuclear enhancement factors of 1.45 and 1.84. The shaded 
areas for the Chamaeleon, R CrA, and Cepheus and Polaris flare spectra indicate the systematic 
uncertainty evaluated in Section [3j 
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Fig. 17.— CO (left) and Av res (right) versus H I emissivities. Each point corresponds to an energy 
bin (see Tables [2] and [3]) . Errors are statistical only. 



Table 4: Masses in the interstellar clouds for each region. 



Region 


/ 




i 




d (pc) 


M co (M ) 


M Avres (M ) 


Chamaeleon 


[295°, 


305°] 


[-20°, 


-12°] 


215 a 


~ 5xl0 3 


~ 2.0xl0 4 


R CrA 


hl°, 


4°] 


[-24°, 


-16°] 


150° 


~ 10 3 


~ 10 3 


Cepheus and Polaris flare 


[100°, 


125°] 


[15°, 


30°] 


300 6 


~ 3.3 x 10 4 


~ 1.3 x 10 4 



Notes." Dame et al. (1987), b Abdo et al. (2010b) 



The obtained X Av values are (0.22 ± 0.01 sta t ± 0.08 sys ) x 10 22 cm" 2 mag" 1 , (0.21 ± 0.01 stat ± 
0.02 sys ) x 10 22 cm" 2 mag" 1 , and (0.14 ± 0.01 stat ± 0.03 sys ) x 10 22 cm" 2 mag" 1 for the Chamaeleon, 
the R CrA, and Cepheus and Polaris flare regions, respectively. With the procedure similar to that 
for CO, we can calculate the mass of additional gas traced by Av Tes expressed as Ma Vi . cs in Table 
H we obtained ~ 2.0xl0 4 M , ~ 10 3 M and ~ 1.3 x 10 4 M & for the Chamaeleon, R CrA, 
and Cepheus and Polaris flare regions, respectively. We thus obtained mass estimates for the 
Chamaeleon and R CrA regions similar to previous ones (Dame et al. 1987) if we consider the total 
mass (traced by Wco and Av res ), although the procedure is not straightforward since the gas traced 
by Aw res is extended in a much larger region of the sky. Detailed study of the matter distribution 
in the interstellar space by comparing 7-rays and other tracers will be reported elsewhere. 



5. Summary and Conclusions 



We have studied the 7-ray emission from the Chamaeleon, R CrA, and Cepheus and Polaris 
flare molecular clouds close to the solar system (< 300 pc) using the first 21 months of Fermi 
LAT data. Thanks to the excellent performance of the LAT, we have obtained unprecedentedly 
high-quality emissivity spectra of the atomic and molecular gas in these regions in the 250 MeV - 
10 GeV range. 

The 7-ray emissivity spectral shapes in three regions agree well with the model for the LIS (a 
model based on local CR measurement), thus indicating a similar spectral distribution of CRs in 
these regions. The emissivities, however, indicate a variation of the CR density of ~ 20 % within 
~ 300 pc around the solar system, even if we consider the systematic uncertainties. We consider 
possible origins of the variation are non-uniform supernova rate and anisotropy of CRs depending 
on the propagation conditions. 

The molecular mass calibration ratio Xqo for the Chamaeleon cloud and the R CrA cloud are 
comparable, whereas that of the Cepheus and Polaris flare region is ~ 2/3 of the others, suggesting 
a variation of Xqo hi the vicinity of the solar system. From the obtained values of Xco> the 
masses of gas traced by Wco hi the Chamaeleon, R CrA, and Cepheus and Polaris flare regions 
are estimated to be ~ 5 x 10 3 Mq, ~ 10 3 Mq, and ~ 3.3 x 10 Mq respectively. Similar amounts 
of gas are inferred to be in the phase not well traced by the Hi or CO lines. Accumulation of 
more 7-ray data, particularly at high energies, and progress in ISM studies, will reveal the CR and 
matter distribution in greater detail. 
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